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ABSTRACT 

In this paper, we revisit the issue of estimating the “fossil” disc mass in the circumprimary 
disc, during the merger of a supermassive black hole binary. As the binary orbital decay speeds 
up due to the emission of gravitational waves, the gas in the circumprimary disc might be 
forced to accrete rapidly and could in principle provide a significant electromagnetic coun¬ 
terpart to the gravitational wave emission. Since the luminosity of such flare is proportional 
to the gaseous mass in the circumprimary disc, estimating such mass accurately is important. 
Previous investigations of this issue have produced contradictory results, with some authors 
estimating super-Eddington flares and large disc mass, while others suggesting that the “fos¬ 
sil” disc mass is very low, even less than a Jupiter mass. Here, we perform simple ID calcu¬ 
lations to show that such very low estimates of the disc mass are an artifact of the specific 
implementation of the tidal torque in ID models. In particular, for moderate mass ratios of 
the binary, the usual formula for the torque used in ID models significantly overestimates the 
width of the gap induced by the secondary and this artificially leads to a very small leftover 
circumprimary disc. Using a modified torque, calibrated to reproduce the correct gap width 
as estimated by 3D models, leads to fossil disc masses of the order of one solar mass. The 
rapid accretion of the whole circumprimary disc would produce peak luminosities of the or¬ 
der of 1-20 times the Eddington luminosity. Even if a significant fraction of the gas escapes 
accretion by flowing out the secondary orbit during the merger (an effect not included in 
our calculations), we would still predict close to Eddington luminosities that might be easily 
detected. 

Key words: accretion, accretion discs - black hole physics - hydrodynamics - gravitational 
waves - galaxies: formation 


1 INTRODUCTION 

The merger of two supermassive black holes (SMBH) is estimated 
to be one of the most intense events of gravitational wave (GW) 
emission and the ability to detect one such event is of paramount 
importance not only to test General Relativity directly but also to 
provide constraints on galaxy formation. 

There is now strong evidence for the presence of SMBHs with 
masses b etween 10^ and 10^ Mp, in the nuclei of the local Uni verse 
galaxies dKormendv & Richstonell9^:lMagorrian et alj 19981: Fer- 
rarese & Ford l2005ll . In a hierarchical scenario for the galaxy evo¬ 
lution, nearby galaxies are the outcome of several mergers between 
progenitors of smaller mass. If each of the progenitors of a galaxy 
merger contains a SMBH in its centre the resulting galaxy will nat¬ 
urally host a pair of SMBHs. 

It is expected that dynamical friction between the binary and 
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the crowded field of stars in which they are embedded can remove 
angular momentum from th e binary and reduce its orbital separa - 
tion down to parsec scales iBegelman. Blandford & Reeslll98d) . 
at which the process is expected to stall du e to the depletion of 
the loss cone jMilosavlievic & M erritfcooi]). Galaxy scale 3D hy- 
drodynamical sim ulations dEscala et al.ll2005l : ICuadra et al.ll200^ 
IPotti et alJuOObh have shown that further shrinkage of the binary 
due to angular momentum loss to the gaseous background can lead 
to orbital separations of the order of 0.1 pc. At separations of the 
order of 0.001 pc, the main driver of orbital decay is the emission of 
GW. The evolution of the system between 0.1 and 0.001 pc is much 
harder to study. 3D and 2D simulations are appropriate to study the 
dynamics of binary-disc interaction on short timescales, and are 
thus well suited to study the final stages of the orbital decay, before 
and during the GW-driv en phase jD’Orazio. Haiman & MacFadyenI 
l2013l:lFatTisetal.l2014li . The overall evolution of the system during 
the long disc-driven dec ay is best followed using simple ID diffu- 
sion models for the disc l lArmitage & Nataraianl2003 : lLodato et alj 
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In this paper, we want to estimate the mass left over in the 
circumprimary disc at decoupling, i.e. when the orbital decay due 
to GW emission becomes dominant. This quantity is the result of 
the long term evolution of the system from large distances (~ 0.1 
pc) down to the decoupling region. For our purposes, thus, it is most 
appropriate to study the problem within a ID diffusion model. 

In previous studies several configurations for the d isc-bi nary 
system have been investigated. lArmitage & NataraianI (l2002h as¬ 
sumed that the binary finds itself embedded i n a gaseous disc tha t 
has already settled down to a steady-state disc. lLodato et alj ( l2009l) 
investigated a configuration in which the binary interacts with a 
finite amount of gas that is brought in the binary vicinity as a con- 
seque nce of the same acc retion episode that has formed the binary 
itself. lOiang et in( l2010l) considered the evolution of a binary em¬ 
bedded in a steady-state disc with a large mass inflow coming from 
the disc exterior. In all these configurations we can describe the sys¬ 
tem as sketched in Figure [T] the accretion disc is rotating around 
a central (or primary) BH of mass Mp, and the secondary BH of 
mass A/s is in Keplerian orbit at a radius a around the primary. 
For simplicity, we assume that the primary BH is located at the 
centre of the accretion disc (thus neglecting the displacement be¬ 
tween the primary BH and the centre of mass of the binary) and 
that the disc is coplanar to the binary orbit ( Ivanov. Papaloizou & 
Polnarev 1l999ll . The disc and the binary exchange angular momen¬ 
tum through tidal interaction in a planet-like dynamics: a secondary 
BH of non negligible mass opens up a gap in the disc at the loca¬ 
tions where the tidal torques equal the viscous torques, with the 
gap following the secondary as it migrates {disc-driven phase). A 
low-mass secondary BH barely perturbs the disc and thus migrates 
on a viscous timescale ty = R^jv, where R is the radial cylindri¬ 
cal coordinate and v is the disc viscosity. Conversely, a secondary 
BH of mass comparable to the disc is able to heavily perturb the 
disc by opening a wide gap but migrates more slowly. Neglecting 
accretion onto the secondary BH, the disc results thus divided into 
two well defined regions: a circumprimary (or inner) disc orbiting 
around A/p and a circumbinary (or outer) disc orbiting around the 

whole binary. _ 

_ The investigations by iLodato et alj (l2009h and 

(201C) lead to apparently conflicting results. Indeed, 

1 20100 estimate that the fossil disc mass is extremely low, of the or¬ 
der o f 10~^ — 10~^ A/(^, de pending on the parameters. On the other 
hand, iLodato et al.l ( l2009h . while not giving specific numbers for 
the inner disc mass, clearly show that the sur face density at decou¬ 
pling can be much higher than estimated bv IChang et alj i2010t) . 
The origin of this discrepancy is unclear, and it could be due to the 
different setup and initial conditions used in the two studies. In this 
paper, we aim to resolve the issue by re-evaluating the inner disc 
mas s, using a disc configu ration close to th e one used bv Chang 
et al. feOlOh . ^ 

As the binary decays further, the angular momentum loss due 
to the emission of GW rapidly increases and eventually takes over 
(GW-driven phase). We define the decoupling time as the time 
when the binary merger timescale becomes shorter than the vis¬ 
cous timescale at the inner edge of the outer disc, i.e. when the 
outer disc is not able to follow viscously the accelerating secondary 
and thus decouples from the evolutio n of the binary (for a detailed 
study on the decoupling process see lFarris et alJlioidh . Assum¬ 
ing for simplicity that the inner disc can drain only by accretion 
onto the primary BH, in the late GW-driven phase the inner disc 
is forced to accrete by the plunging secondary BH that migrates 


Qian^_et_^ 

Chang et al. 



Figure 1. Sketch of the binary system. The accretion disc, that orbits around 
the primary BH of mass Mp is divided into two regions by the presence of 
the secondary BH of mass Ms that orbits in Keplerian motion at a radius a\ 
a circumprimary disc (inner disc) and a circumbinary disc (outer disc). The 
inner disc is truncated at Sedge whose size depends on the intensity of the 
tidal torques. 


much more rapidly. The extremely rapid accretion of the whole in¬ 
ner disc (whose mass is substantially frozen during the GW-driven 
phase) might lead to a suner;Eddington_flarej_^oteJhaLBaruteau, 
Ramirez-Ruiz & Masset T2012l) . using 2D simulations, show that 
such rapid accretion can be reduced due to the funneling of the 
inner gas out of the secondary orbit during the GW-driven phase, 
estimating that up to 80 percent of the inner disc mass can escape 
accretion in this way. Such results, however, need to be confirmed 
in 3D simulations and using a w ider range of parameters, especially 
for th e disc aspect ratio (which iBamteau. Ramirez-Ruiz & Massed 
I2OI2I fix to a relatively large value). 

The paper is organized as follows. In Section |2] we illustrate 
the physical model that we adopt to describe the disc-binary system 
with a detailed discussion of the torque implementation. In Sec- 
tion[^we describe the results of a set of simulations with different 
parameters and the overall evolution of a fiducial run. In Section|4] 
we discuss our results in comparison with those in literature and 
we propose an analytic argument to explain the discrepancy in the 
estimate of the inner disc mass. In Section |5] we present our con¬ 
clusions and the outlook for this work. 


2 PHYSICAL MODEL 

In order to describe the coupled evolution of the disc-binary system, 
we have implemented a time-dependent ID model that solves for 
the disc surface density S(R,t) and for the binary separation a(t) in 
a self-consistent way. 


2.1 Disc dynamics 


The evolution of the disc surface density in presence of the sec¬ 
ondary BH is descri bed by the classical hydrodynamical equation s 
of accretion theory dLvnden-Bell & Pringld l l97^ |Pringlg[l98lh : 
the continuity equation 


dY. Id 


( 1 ) 























































3 


and the angular momentum conservation 


I ^ ^ ^ ^ + AtE, (2) 

where v = (vr, V(p,vz) is the fluid velocity, Gy = 2kvLR^Q! is the 
usual viscous torque and Ax is the specific tidal torque exerted by 
the binary on the disc. Using (D and © we obtain an expression 
for vr and substituting it back in © we obtain the time evolution 
equation: 


dL 


3 d 

rJr 




2 d [AjE 

rJr 


(3) 


where v{R) is the disc viscosity and Q.{R) = ^GM^/R? is the disc 

angular velocity. _ 

Making use of the impulse a pproximation dUin & Papaloizoul 
ll979al : [Binnev & Tremaineiri987n we can compute a simplified ex¬ 
pression f or At(/^ 1 that has been widely used in literature (Lin & 
Papaloizou 1 1979bl: Armitage & NataraianI 20021 : Lodato & Clarkel 
|2004 iLodato et al jboOOll: 


/ 


Aj = ^Mr^[-- 


At = n 


R <a, 


R> a, 


(4) 


(5) 


where / is a normalization factor and A = R~a. The change of 
sign across R = a in the tidal torque accounts for the behaviour 
of the secondary BH that removes angular momentum from the 
gas inside its orbit and adds angular momentum to the gas outside. 
The net effect is that the satellite repels the gas from its orbit, thus 
creating an annular gap whose width depends on the intensity of 
the tidal torques with respect to the pressure and viscosity gradients 
that act to close the gap. 

Since theaboveex£tessira^divergesJn_g^^_aj^llgwmg_Lin & 
Papaloizou (Il986h and lSver & Clarkel d 19951) we smooth the torque 
by taking: 


A = max{//,7?Hj i 


( 6 ) 


where H is the disc thickness and = a(g/3)*^^ is the Hill’s 
sphere radius of the secondary BH. This smoothing prescription 
follows from the basic idea that the gap should be larger than both 
H (otherwise it would be replenished by the pressure gradients in 
less than a dynamical timescale) and R^ (that should be depleted 
due to the accretion on the secondary bhIII 

It is easy to see that the intensity of the tidal torque determines 
the gap width since the gap edges are the locations where the vis¬ 
cous torques, that try to refill the gap with a diffusive behaviour, are 
balanced by the tidal torques that act to keep the gap open. Analo¬ 
gously to Ax, we can define the specific viscous torque Ay dividing 
the annular viscous torque Gy by the mass of the annulus, namely: 


Gy 

IkRILA 



3 SIR 
2~a'' 


(7) 


where the last equality holds for a Keplerian disc. The gap width 
is thus the value of A for which Ay = Ax. In Figure 2 we illustrate 
how, if the intensity of Ax changes, then the location of the gap 
edges changes as well: for a given viscous torque intensity (that in 


* We should note that apaid from this smoothing prescription we have not 
taken into account the accretion o n the secondary B H in the model. This is 
a good approximation for ^ = 0.1 jparris et alj2014h . hut some accretion is 
expected for q > 0.3, which might slow down migration somewhat. 



Figure 2. Different torque intensity induce different gap sizes. Two torques 
of different intensity are shown, the upper one being more intense than the 
lower one. The gap inner edge (Sedge) is located where the tidal torque Ax 
equals the viscous torques Ay (that here for simplicity has been fixed to be 
constant). It can be seen that, for a given Ay, the gap edge produced by 
the stronger torque is located at R"^g^, that is smaller (i.e. the gap is wider) 
than that produced by the weaker torque (S^^g^). The same result applies 
for R> a. 


the figure is assumed constant for simplicity) a stronger tidal torque 
truncates the inner disc at a smaller radius. Therefore, in order to 
implement a realistic tidal torque, the parameters describing the 
torque (/ and the smoothing prescription) have to be chosen care¬ 
fully, so as to reproduce the gap size obtained by more sophisticated 
3D simulations. 

Since the size of the gap in the disc, as well as the migration 
rate of the secondary, depends on the intensity and the spatial distri¬ 
bution of the tidal torques, a correct smoothing of the tidal torques 
is necessary not only to avoid numerical issues, hut mostly to avoid 
an unphysical estimate of the interaction between the disc and the 
binary that can lead to a wrong estimate of the inner disc mass. 


2.2 Tidal torque smoothing 

As shown bv iGoldreich & Tremaind i I 1979L[1983) the disc and the 
secondary interact hy exchanging angular momentum at discrete 
resonant locations, the so-called Lindblad Resonances, where a 
characteristic frequency of the secondary matches an epicyclic fre¬ 
quency within the disc: at these locations the disc responds to the 
presence of the secondary through the excitation of density waves 
that carry energy and angular momentum in the disc and dissipate 
through shocks. 

In a Keplerian disc the Lindhlad Resonances are located at 

radii: 

/ 1^/3 

RL = flfl±-j ra=l,2,... (8) 

which shows that for low m the Lindhlad resonances are 
far from the position of the secondary whereas for increas¬ 
ing m they accumulate towards R = a. We find that the in¬ 
nermost Lindblad Resonance (hereafter, IMLR) is placed at 
Rimlr = ~ 0.63a, whereas the outermost one (OMLR) is 

placed at Rqmlr = 2^/^a ~ 1.59a. Since the IMLR and the OMLR 
are the farthermost resonances from the satellite, they constitute the 
boundaries of the region that can effectively contribute to the tidal 
torque: outside this region, the disc can exert no torque on the satel- 
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Table 1. Inner and outer gap edges obtained by lArtvmowicz & Lubowl 
h994 performing 3D SPH simulations, for « — 0. We report the widths 
Mmlr and v'^OMLR of the gaussian smoothing that is needed to reproduce 
the same gap edges. 



Gap Edges 

Gaussian widths 


Inner 

Outer 

vriMLR 

Womlr 

0.11 

0.46a 

1.8a 

370// 

75// 

0.43 

0.38a 

1.8a 

370// 

83// 


lite, nor the satellite can exert any torque on the disc. As noted in 
iLodato et alJ bOOgI) . since At scales as q^, for high q it can oc¬ 
cur that the analytical expressions and 0 gives a non negligi¬ 
ble contribution also outside the region of LRs. In order to avoid 
this unphysical contribution, we need to smooth out the torque 
for R < Aimlr and R > Aqmlr- We can do this by fine-tuning the 
torque formulae 0 and 0 in two different ways, depending on the 
mass ratio of the binary. 

For ^ <C 1 the torque formulae 0 and 0 already gives no 
contribution outside the IMLR and the OMLR, and the gap size is 
determined by the parameter /. Therefore, in this case, it is suffi¬ 
cient to fi ne-tune the value of f to reprod uce realistic gap sizes. 
Following lArmitage & NataraianI | |2002[) that performed high- 
resolution 2D simulatio ns using the ZEUS hydrodynamic code 
dStone & Normanlll992h . we take / = 10 ^ with which we repro¬ 
duce gaps of approximately the correct size. 

For larger mass ratios, the torque is still significant outside the 
LRs region and thus we smooth it out with a gaussian cutoff by 
setting: 


Ai = -^q^d}R^ 


/rY 

/ /?-/?imlrV 



va; 

V Wimlr / 


R ^ ^IMLR 


(9) 


At = 


-q^d}R^ 



f R — f^OMLR 

V ’^OMLR 


2' 


R ^ f^OMLR 


( 10 ) 


where wimlr and vvqmlr determine the sharpness of the cutoff 
and for Rimlr <R< Romlr the formulae 0 and 0 still hold. 
We fine-tune the values of wimlr and womlr so that the gap 
edges (measured after the secondary has opened a clear gap in the 
disc) are located in the same positio ns computed by Artymowicz 
& Lub ow (Il994h with .3D SPH simulations. lArtvmowicz & Lubowl 
(I1994 studied circumbinary discs around binaries with g’ = 0.11 
and q = 0.43, on circular and eccentric or bits, and for discs with 
differe nt Re ynolds numbers. Furth e rmore , lArtvmowicz & Lubowl 
(Il994h and fPapaloizou & Pringld (Il977h showed that for zero- 
eccentricity binaries the location of the gap edges does not depend 
on the value of the Shakura-Sunyaev viscosity parameter a (see Eq. 
Ell. In Table [T] we report the values of wimlr and vvqmlr that we 
obtain for each mass ratio. Since the binary we are considering for 
our purposes is on a circular orbit and the values found for wimlr 
and Womlr weakly depend on the mass ratio, for our simulations 
we always adopt the values wimlr = 370/7 and wqmlr = 75//. 

Note that the formula used here for the local torque is clearly 
an approximation, used in our ID model to reproduce the correct 
gap width^ndsecond^^jni^rationj^ej^egentji^toionsb^Far- 
ris et al. i | 2014 allow to compute the local torque as a function of 
radius. The torque obtained in this case has a different functional 
form from the one used here. In particular, it has an oscillatory be¬ 


haviour and does not vanish until 2-3a. However, for comparable 
values of q, Farris et al co nfirm the gap widths found by Artymow¬ 
icz & Lubow jl994- The exact shape of the torque close to the 
secondary does not affect any of our results (indeed, in many cases 
one can treat the secondary, simply imp osing an effective bound¬ 
ary condition, see lSver & Clarkal 19951) . as long as the gap width 
and migration rate are reproduced correctly. If we had kept our an¬ 
alytic al torque formula, allowing it to extend to 2-3a jFarris et alj 
|2014 . we would have obtained an incorrect gap width. Hence we 
have preferred to be consistent with our approach and smooth the 
torque in order to reproduce reasonable gap widths. 


2.3 Disc energetics 

In an accretion disc surrounding a SMBH binary we expect to find 
regions characterized by different regimes, depending on the dis¬ 
tance from the central BH, the tem perature, the d i sc thic kness, etc. 
Following the approach adopted in lLodato et ^ ^20091) . our code 
is capable of reproducing both radiation-pressure and gas-pressure 
dominated regimes (we expect the first to occur in the inner part of 
the disc and the latter in the outer ones), and also to take into ac¬ 
count both Thomson scattering and free-free absorption as sources 
of opacity. 

In order to compute the viscosity of the disc, we adopt the 
standard a—prescription jShakura & Sunvae^ll973h : 

v = ac,H, (II) 


where Cj is the sound speed and a is a dimensionless number 
smaller than unity. In ord er to avoid viscous and therm al instabili¬ 
ties in the disc, following iLightman & EardlevI (Il974 we assume 
that the viscous stresses are proportional to the gas pressure (and 
not to the total pressure, sum of gas and radiation contributions), 
and therefore Cj = with being the Boltzmann con¬ 

stant, Tc the midplane temperature, /i = 0.67 the mean molecular 
weight (ionized Hydrogen) and nip the proton mass. 

The midplane temperature profile is computed by solving the 
vertical energy balance at each radius. First, we suppose that the 
energy dissipated by the viscous processes (released mainly in the 
disc midplane) is completely irradiated by the disc at its surface by 
imposing: 

^vzn2 = (7sBr/, (12) 

o 

where (Tsb is the Stefan-Boltzmann constant and Ts is the surface 
temperature. Then, we assume that the disc is optically thick to its 
own radiation, i.e. the optical depth t = k'Z/ 2 3> 1, where K is the 
opacity (given by the dominant source among Thomson scattering 
or free-free absorption). Finally, plugging the a—prescription O 
in Eq. d and using the fact that = 47’c/3'r'/^, we find that the 
midplane temperature is given by: 


Tc = 


/ 9 akB 
V 16 RmpCTsB 



1/3 


(13) 


Once computed the viscosity and the central temperature, the disc 
thickness can be calculated from the vertical hydrostatic equilib¬ 
rium that results from the balance between gravity and the total 
pressure. At each radius the structure of the disc is fully determined 
by V, Tc and //, and we use an iterative procedure that converges to 
consistent valu es for these three par ameters. 

Note that iLodato et alJ l l2009li have shown that, in order to 
simultaneously conserve energy and angular momentum in the ID 
formulation that we adopt here, a tidal heating contribution should 
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be added at the disc edge to account for the dissipation of non- 
axisymmetric structures induced by the tidal force. W e neglect here 
this co ntrihution, mostly to compare more directly to lChang et alJ 
( l2^ (who also do not include it). This tidal heating term could 
mo dify some of the results in limite d portion of the parameter space 
(see lKocsis. Haiman & Loebl2012h . when the gap is not fully open. 


2.4 Binary dynamics 


To derive the equation for the binary evolution we simply impose 
that the angular momentum of the whole disc-hinary system is con¬ 
served and we take into account the angular momentum lost due to 
the emission of GW. The binary separation thus evolves as follows: 


4n 


r 

J Rin 


AjLRdR 


8c 



(14) 


where f2s = y GM-p/a^ is the angular velocity of the secondary and 

Rs is the Schwarzschild radius of the primary, R^ = IGMpjc^. The 
first term on the right hand side is the tidal back reaction of the 
disc on the binary and the second term is due to the GW emission 
dPeters & Mathew3ll963HPetersll9^ . 

We note that while the direction of migration due to the tidal 
torque can change in time since it depends hoth on the intensity of 
the torque and on the surface density distrihution, the emission of 
GW always acts to shrink the binary. 


3 RESULTS 


3.1 Numerical code: setup and initial configuration 

We solve numerically the coupled equations ([^ and ( 114b with a 
ID finite-difference numerical scheme that allow us to compute the 
time evolution time of and a{t). We use a radial logarithmic 

grid of 200 points, with the inner radius equal to the innermost sta¬ 
ble circular orbit (ISCO) around the central BH, i.e. = 3Rs, and 
the outer radius Rout = lOao where oq = lO^i^s is the initial posi¬ 
tion of the secondary BH. At we apply a zero-torque boundary 
condition to implement the accretion on the central BH, whereas at 
^out we can choose to set either a zero-torque boundary condition 
or an inflow of material Mext. depending on the scenario we want 
to model. _ _ 

In previous investigations, lArmitage & NataraianI l l2002h 
initialized the disc with a steady-state disc corresponding to 
Mext = 13T0yr^*, while during the simulation they assumed 
iHext = 0. With this choice, the secondary BH results parachuted 
inside a disc that has already settled down to a steady state and 
therefore the binary hardening is driven by the almost infinite reser¬ 
voir of gas upstrea m and not by an inflow of material coming from 
Rout- More recently. lUodato et al.l l l2009t) investigated a scenario in 
which the SMBH binary interacts with a finite amount of gas that 
evolves together with the binary and is far from a steady state; as 
initial configuration they assumed a constant surface density pro¬ 
file between 0.8ao and 2ao and Mext = 0 to simply model this sce¬ 
nario. We performed several tests of our code using these two con¬ 
figurations andjveMe_able_^o_re 2 roduce_the^esults_of^Armitage 
& Nataraian j2002[l and Lodato et alj ( l2009h to a high degree of 
accuracy. 

Despite being less realisti c than the previous s cenario, to com¬ 
pare our results to the work of IChang et alj ( l2010h . in all our sim¬ 
ulations we choose initialize the disc to a steady-state disc of low 


mass. To do this, instead of adopting analytic solutions, we have ob¬ 
tained the steady state configurations by solving numerically with 
our code Eq. © without the tidal term on the rhs and applying 
a Mext 7 ^ 0. To probe the accuracy of our code, we successfully 
proved that the surface density profile of this initial steady state 
is in good agreement with the standard solutions of the Shakura- 
Sunyaev a—disc. The initial accretion rate through the disc is thus 
constant and is set to a very low value, much smaller that the Mext 
that we then apply. As the simulation starts, the disc starts to fill up 
with mass and the accretion rate thus starts to slowly grow, as the 
disc tries to evolve toward the new configuration appropriate for the 
higher accretion rate. At the initial time-step, we set the position of 
the secondary BH in R = ao and we let the system evolve: at the 
beginning, the secondary BH clears a gap in the disc and, once the 
tidal interaction becomes strong enough, it starts migrating. 


3.2 Resulting fossil disc mass 


How much mass is left in t he inner disc at decoupl i ng ha s been the 
subject of debate recently. lArmitage & NataraianI (|2002[) estimate 
a relatively large m ass based on very idealized initial conditions. 
lUodato et alJ ( l2009h . using more realistic conditions found a rel¬ 
atively large mass, too, with super-Eddington accretion rates once 
the inner disc is rapidly accreted du ring the GW-driven binary de¬ 
cay. Conversely. IChang^ 'Zl dTOld) found that the fossil disc mass 
is extremely low, being of the order of planetary masses. Thus, our 
main aim is to resolve the discrepancy and estimate what is the 
m ass of the fossil disc at decoupling, comparing our results to those 
of IChang et al.l ( l2010l) . We thus p erform a number of simulations, 
using the very same para meters of Chang et al.| ( l2010l) . Our code is 
similar to the one used bv ichang et al, 1 2010l) in many respects and 
we would thus expect to reproduce their results. 

In the simulations that we perform we are interested in com¬ 
puting the evolution of the accretion luminosity Lacc as the bi¬ 
nary approaches to the merger and we do that by simply com¬ 
puting Lacc = rjMc^ where T) = 0.1 is the accretion efficiency and 
M is the instantaneous mass accretion rate onto the primary BH 
(measured as the mass per unit time flowing inside the ISCO, i.e. 
R < Rin). As described above, we call decoupling radius the loca¬ 
tion where the merger time driven by the emission of GW becomes 
smaller than the merger time due to the tidal interaction with the 
disc. At that point, the binary decouples from the dynamics of the 
disc and evolves on a ti mescale that is too short for the disc to re¬ 
spond viscously (but see lEarris et al.ll201^ . During this phase, the 
mass of the inner disc is substantially frozen since the secondary 
forces it to accrete onto the primary at a rate much higher than 
the viscous one. During the evolution of the system, we compute 
Lacc from M at each timestep and we estimate the peak luminosity 
Lpeak = max(Lacc) that is reached immediately before the merger. 
We can put a lower limit to Lpeak estimating the average luminos¬ 
ity that is produced by the accretion of the inner disc in the time 
interval between the decoupling and the merger, namely 


-'peak 




Md,in (tdec ) 
^merger idee 


rjc 


(15) 


where M^^n (f) is the time-dependent mass of the inner disc, defined 
as follows: 


Md,in(r)= r°''°27tRI.{R,t)dR, 

J Rin 


(16) 


where i?edge is the outer edge of the inner disc, as shown in Eigure[T] 
Clearly, this estimate is valid if we assume that the whole inner disc 
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is acc reted during the merger tearuteau. Ramirez-Ruiz & Massed 

I2OI2I) . 

We perform a series of simulations exploring differ¬ 
ent regions of the parameters space by varying the primary 
BH mass Afp = 10^, 10^, 10*Af©, the external accretion rate 
A?ext = 0.1, O.OlAfgyr^*, the mass ratio ^ = 0.3, 0.1 and the 
Shakura-Sunyaev viscosity parameter a = 0.1, 0.01, whereas we 
keep fixed the initial separation aq = lO^i^s and the initial disc 
mass A/(j = 10* M©. 

The results of the simulations are shown in Table We re¬ 
port the inner disc mass at decoupling and the peak value of 
the accretion luminosity Lpeak just prior to merger. In almost all the 
cases we obtain peak luminosities that clearly exceed the Edding¬ 
ton luminosity and in four cases the flare is super-Eddington by one 
order of magnitude. It is remarkable that even though A/j is rel¬ 
atively small (in many cases a few percent of one solar mass), its 
forced accretion during the final inspiral of the secondary results in 
an extremely luminous flare. 

The last two columns in Table 2 show t he fos sil disc mass and 
peak luminosity as obtained bv ichang et alj l l201(]ll . There is a strik¬ 
ing difference in the estimated disc mass by up to three orders of 
magnitude with respect to ours. This systematic difference is found 
for all parameters used and results also in a significant difference in 
peak luminosity. We discuss the origin of this discrepancy in sec- 
tion|4l after describi ng the dynamical ev olution of a fiducial run. 

Note that while lchang et Ii]( l2010l) obtain an almost linear re¬ 
lation between the fossil disc mass and the external accretion 
rate A?ext, in our case the dependence on Mext is much weaker. For 
most of our simulations (with the exceptions of the low viscosity 
and high primary mass cases), the inner disc mass scales as A?ext • 
Increasing A/ext leads to a faster migration rate which implies that 
the inner disc has less time to drain out, but a linear scaling between 
A/(] in and Aifext is not necessarily expected. Indeed, the dependence 
of the migration rate on the external accretion rate is non linear as 
it depends both on the surface density and the temperature at the 
inner edge of the outer disc, thus making it not straightforward to 
find a simple prescription to relate these two quantities. 

3.3 Evolution of a fiducial run 

Let us describe the evolution of the system for a fiducial case with 
A/p = 10* A/©, g = 0.1, A/ext = 0.1 A/© yr^^ and a = 0.1. 

We choose to initialize the disc with the same choice of Chang 
et al. 1 I2OIOI) . i.e. with a low mass disc of mass 10* A/,^. appropri- 
ate for a Shakura-Sunyaev disc accreting at a very low rate, e.g. 
10^® — 10^*A/©/yr. Note that the purpose of this initial config¬ 
uration of the disc is to give a simple realization of the planet¬ 
like regime in which the binary evolu tion takes place an d it can 
be shown (for a detailed discussion see lChang et akll^OlOl) that the 
inner disc mass at decoupling - and thus the peak luminosity just 
prior to merger - does not depend on the initial mass of the disc. 

The ge neral evolution of th e system is very similar to the one 
described in ichang et alj j2010l) . Following Figure[3 in each snap¬ 
shot we show the surface density profile (solid line), the viscous 
timescale (dashed line) and the merger timescale /merger = a/a (dot¬ 
ted line). Time evolution proceeds as follows: top-left, top-right, 
bottom-left, bottom-right. At the beginning (Figure[3] top-left), the 
system reaches an initial quasi-stationary state in which the outer 
disc accumulates some mass due to the external accretion rate but 
the tidal torques are still not strong enough to push the secondary. 

Starting from this moment, the outer disc increased surface 
density begins to influence the secondary and to drive its inward mi- 


Table 2. Comparison of simulations results. We perform simulations 
for different values of Mp = 10^, 10*, 10* A/©, Mext = 0.1, 0.01 M© yr“’, 
^ = 0.1, 0.3 and a = 0.1, 0.01. In the Table, masses are in units of M© and 
luminosities in units of LEdd- We report the inner disc mass at decoupling 
A/dpn and the luminosity peak just prior to merger . Fo r better compai'- 
ison, we report the values obtained bv IChang et al.l I2010l) in the last two 
columns, and d*\. 

’ cl,in peak 


Mp 

1 

^ext 

a 

^d,in 

^peak 


^peak 

10® 




xlO^* 


xlO^* 



0.1 

0.1 

0.1 

0.95 

1.97 

0.0012 

0.146 




0.01 

17 

13.7 

0.012 

0.277 



0.01 

0.1 

0.41 

0.86 

0.00012 

0.015 




0.01 

8.4 

6.44 

0.0012 

0.029 


0.3 

0.1 

0.1 

1.47 

1.7 

0.0001 

0.038 




0.01 

25 

8.8 

0.001 

0.023 



0.01 

0.1 

0.68 

0.92 

0.000008 

0.004 




0.01 

13 

5.32 

0.0001 

0.002 

10* 




xlO^* 


xlO^* 



0.1 

0.1 

0.1 

183 

6.22 

0.66 

0.154 




0.01 

2970 

20.8 

46 

1.077 



0.01 

0.1 

83 

1.38 

0.073 

0.018 




0.01 

1600 

9.66 

0.64 

0.154 


0.3 

0.1 

0.1 

300 

2.9 

0.062 

0.015 




0.01 

450 

23.5 

0.56 

0.115 



0.01 

0.1 

140 

1.51 

0.0062 

0.001 




0.01 

2400 

8.22 

0.075 

0.015 

10* 

0.1 

0.1 

0.1 

19 

2.67 

0.2 

0.5 




0.01 

34 

0.47 

0.9 

2.154 



0.01 

0.1 

11 

1.31 

0.029 

0.069 




0.01 

35 

0.48 

0.18 

0.438 


0.3 

0.1 

0.1 

26 

2.15 

0.024 

0.052 




0.01 

34 

0.189 

0.15 

0.323 



0.01 

0.1 

22 

1.24 

0.0033 

0.007 




0.01 

34 

0.19 

0.033 

0.069 


gration. The gap edges are the locations where the viscous torques 
equal the tidal torques and, since the latter depend on the position 
of the secondary, the gap edges move together with the secondary. 
As a result the system self-regulates in such a way that the vis¬ 
cous timescale at the gap edges equals th e merger timescale during 
the disc-driven phase dChang et al.ll201ol) . Then, as the disc-driven 
migration continues (Figure 0 top-right), the gap moves inwards 
along with the secondary, as we expect for a Type II migration. 
The orbital radius of the secondary has decreased by one order of 
magnitude, and the inner disc has a slightly enhanced surface den¬ 
sity. The viscous timescale at the gap edges continues to match the 
merger timescale (Figurebottom-left). When the binary separa¬ 
tion has decreased by another order of magnitude (a ~ 100/?s) the 
emission of GW rapidly takes over the tidal interaction and starts 
driving the migration. At this stage the outer disc is substantially 
decoupled from the evolution of the binary and the inner disc is 
pushed inwards more and more rapidly as the binary shrinks. Even¬ 
tually (Figure 1^ bottom-right) the increasing velocity of the sec¬ 
ondary (accelerated by the growing GW emission) forces the inner 
disc to accrete on a timescale much shorter than the viscous one, 
thus producing the spike in the surface density profile and the con¬ 
sequent flare in luminosity. In the last snapshot /merger is not shown 
because it is much smaller (about days) than the viscous one. 

In Figure |4] we show the accretion luminosity on the primary 
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Figure 3. Disc and binary evolution for the fiducial case with Mp = 10^ Mq, q = Q.l, =0.1 Mq yr“* and a = 0.1. Each snapshot shows the surface density 
profiles 21 (solid lines), the viscous timescale profile (dashed lines, matched to the right axis), the merger time (merge (horizontal dotted line, matched to the 
right axis) and the secondary position (marked with a black dot). The snapshots have been taken respectively at 5 Myr (top-left), 12 Myr (top-right), 12.941 
Myr (bottom-left) and 12.948 Myr (bottom-right) after the beginning of simulation. 


BH as a function of time during the overall evolution of the sim¬ 
ulation. During the first 10^ years the accretion rate is very low, 
and it reflects the initial condition for the disc surface density (see 
Section [3}. During this time the disc mass huilds up as a conse¬ 
quence of the imposed external accretion rate, slowly increasing the 
disc luminosity to values closer to typical AGN values. As the disc 
mass grows, it pushes the secondary inwards, until, after roughly 
10^ years, the secondary migration becomes dominated hy grav¬ 
itational wave emission. Finally, it can he clearly seen the super- 


Eddington flare due to the extremely rapid accretion of the inner 
disc that occurs at the end of the GW-driven inspiral. 


4 DISCUSSION 

4.1 Effects of tidal torque smoothing on luminosity peak 

Comparing our results with those hy IChang et al.l 1 I 2 OIOI) we see 
that these authors predict much lower inner disc masses and, conse¬ 
quently, much smaller peak luminosities for all the computed mod- 
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t/yr 


Figure 4. Accretion luminosity on the primary BH as a function of time 
during the overall evolution of the disc-binary system for a fiducial run 
with Mp = 10^M 0 , q = 0.1, Mcxt, 01 = 0.1. The SMBHs coalesce around 
13 Myr, immediately after the apparent super-Eddington flare due to the 
rapid accretion of the inner disc. 


els by up to three orders of magnitude. For a better comparison, we 
report their results in the last two columns of Table [21 In order to 
explain the discrepancy in the estimate of the inner disc mass at the 
decoupling, we have performed several tests before finding that the 
issue was in the implementation of the torque. 

Before discussing the torque implementation, we would men¬ 
tion that a possible source of discrepancy may arise from the differ¬ 
ent viscosity calculation: as explained in Section [23l while we ne¬ 
glect radiation pressure when computing viscosity in order to main¬ 
tain the disc stable, we take it into account when computing the disc 
thickness assuming that the disc is support ed by both gas and radi- 
ation pressure. Therefore, differently from IChang et alj bOlOl) that 
completely neglect radiation pressure, our code is able to describe 
radiation pressure dominated regions as well. During our tests, we 
see that the radiation pressure actually becomes dominant only in 
the inner disc in the very late phases just prior to merger and in the 
inner edge of the outer disc after a significant amount of gas has 
piled up. However, the estimate of the inner disc mass seems not to 
be affected by this choice. 

Another source of discrepancy might come from the way the 
gap is initially opened in the disc: while our code computes the gap 
opening consistentl y with the torque calculated at each timestep, 
IChang et alj ( 1201 Oh artificially open a cavity in the disc between 
O.SflQ ^ ^5 2ao in order to simply model the gap clearing. At the 

very beginning of the simulation these different implementations 
produce visible differences in the gap width, however after a few 
Myr the gap width results determined by the intensity and the spa¬ 
tial distribution of the tidal torque rather than by the initial clearing. 

We thus conclude that the source of the discrepancy is in the 
implementation of the torque. Indeed, while for small mass ratios 
the LRs responsible for the gap formation are those closer to the 
satellite (and therefore the torque intensity is determined by the 
parameter /), for the larger mass ratios q^OA that we are using 


here, the torque given by Eq. Hi and ([5} is a significant overesti¬ 
mate, and thus needs to be smoothed out as we described in section 
12.21 above. In Fig. [5] we show the evolution of our fiducial model 
with Mp = \QpM q, q = 0.1, Mext = O.lMQyr^* and a = 0.1, but 
now computed without smoothing the torque outside the Lindblad 
resonances. Note th at the evolution sho wn here is almost identical 
to the one shown bv ichang et al.l ( I2OIOI1 . 

Comparing the snapshots in Figure|5]with those in Figure|3]it 
is apparent that, without smoothing the torque, the code produces 
a larger gap (the width being approximately 0.94a instead of 0.7a) 
and smoother gap edges. Such wide gaps are an artifact of the ID 
model with no smoothi ng and are not consisten t with the more de¬ 
tailed 3D simulations bv lArtvmowicz & Lubowl ( Il994h . As a result, 
starting from the initial stages (Figure [S] top-left) the inner disc is 
truncated at a much smaller radius (as we expect after the discus¬ 
sion in Section lT2l . As described in ichang et alj(l2010h the surface 
density of the disc regulates in such a way to keep the viscous time 
at the outer edge equal to the migration tim^^. If the outer edge 
of the inner disc is located much further in, this requires a much 
smaller surface density, thus leading to a significant underestimate 
of the inner disc mass. 

In this way, we can provide an analytic argument to quantify 
the extent to which the inner disc mass is underestimated if the 
torque is not smoothed. At the beginning of the simulation, when 
the binary separation is large (a ~ lO^Rs)! in ths region where the 
inner disc is truncated by the tidal torques exerted by the satellite 
(at R ~ Redoe)^ ths disc is gas pressure dominated and opacity is 
dominated by Thomson scattering, and thus we can write viscosity 
as 


V = 


\^ilnpJ yiecMpassJ 


(17) 


where we have substituted Tc from Eq. l ll3t in the a—prescription 
([IB and used H = Cs/D that holds for a geometrically thin disc. 
With this expression for viscosity, we can estimate the viscous 
timescale at the inner gap edge: 


^v(7^edge) ' 


p2 

^edge 


^edge 


V(7?edge) £2/^ ■ 


(18) 


We have seen that during the secondary migration the gap moves 
together with secondary: this means that the viscous time at /?edge 
has to match the secondary migration time (Fig. [5ji ,b,c). From the 
above equation we see that, if the inner disc is truncated at a smaller 

radius, then the surface density has to adjust itself to a lower value 

3 /2 

^(■^edge) ^edge’ 1° 8*''® tv{Redge) = fmerger- To a first approxima- 
tion we can then compute the inner disc mass: 


^d.in 



iTtRLdR ~ , 


(19) 


where the strong dependence on i?edge reflects the dependence of 
viscosity on surface density in the gas-pressure dominated region. 
From the above result we can now estimate the discrepancy in the 


^ Note that iDuffell et alj 1201 4h and iDurmann & KlevI 1201 4|) have per¬ 
formed numerical simulations of Type II migration where the migration 
timescale did not exactly match the viscous timescale in the disc. This might 
be an interesting phenomenon but needs to be studied with greater detail, 
and in particular it is not clear whether such migration rates can be sup¬ 
ported over long timescales, where a global balance of angular momentum 
could be established. A full discussion of this issue is clearly beyond the 
scope of the present paper. 
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Figure 5. Disc properties computed without torque truncation for the fiducial case with Mp = q = 0.1, ilifext = 0.1 M@ yr“* and a = 0.1. Each snapshot 

shows the surface density profiles Z (solid lines), the viscous timescale profile ty (dashed lines, matched to the right axis), the merger time fmerge (horizontal 
dotted line, matched to the right axis) and the secondary position (marked with a black dot). The snapshots are taken at the same binary separations of Figure[3] 
Time evolution proceeds as follows: top-left, top-right, bottom-left, bottom-right. 


inner disc mass computed by two different torque implementations. 
Adopting the same notation as in Figure |2l and are the 
locations of the gap edges produced respectively by our code (with 
torque smo othing outside the L Rs) and by the torque implementa¬ 
tion used in ichang et alJ OOlOt) (without torque smoothing). Flence 
we obtain that the ratio between the inner disc masses computed in 
the two cases is: 


For the fiducial case that we are considering, we get (at the time 
at which a gap has fully developed but the secondary has not yet 
started migrating) ^ 4.66 x 10^rg and ~ 1.12 x 10^rg 
thus leading to: 


"d.i„ 

Kin 


150, 


( 21 ) 


K.m pedgeV^' 
Km " 'v^dgej 


comparable to the discrepancy that we found between our results 
and those in IChang et al.l 1 I 2 OIOI 1 . We have verified that this ratio is 
not strong function of time. 
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4.2 Limitations of the present calculations 


5 CONCLUSIONS AND OUTLOOK 


Some remarks should be made about the assumptions of the code, 
and the validity of the results. In our code the inner disc can only 
accrete on the primary BH, forced by the GW-driven inspiral of the 
secondary, thus neglecting any possible mechanism that is able to 
drain the inner disk. Therefore, the estimates of the peak luminosity 
that we report should be regarded as upper limits that hold in the 
optimistic case in which there is no inner disc leakage, e.g. across 
the gap or through jets. 

Some processes that we are not considering in our simulations 
might be responsible for a reduction of the inner disc mass in the 
late phases of the binary shrinking. Indeed, given the ID nature 
of our code, by definition we are not able to describe off-plane or 
not cylindrically-symmetric gas flows. However these flows, that 
would require 2D or, better, 3D SPH simulations to be investigated, 
might become relevant only in the innermost regions of the disc, 
where the tidal interaction and the accretion on the primary cause a 
significant heating of the gas, possibly allowing the disc to become 
geometrically thick. 

Another possible leakage of the inner disc may be due to a 


Masset |2012 

). Gas flow across the san has been studied bv Brvden 

et al. ill999h; 

Artvmowicz & Lubo^|l996|) and. more recently, bv 


Fadven HOlsUFarris et al.ll201^ . In general, since the intensity of 
the torque scales with q^, we expect that a low-mass secondary 
is less efficient in producing a clear gap, thus allowing a certain 
amount of gas to be funneled through the gap. As the secondary 
mass increases, the mass flow through the gap decreases signifi¬ 
cantly, reaching a minimum value of the order of 10 percent of 
the unperturbed_accrefion_^ate_Jor_ 2 _^£J^D’Orazio, Haiman & 
MacFadven l2013li . For even lareer secondary masses, the cavity is 
significantly deformed, becoming elliptical and allowing gas to en¬ 
ter the cavity. Thus the mass flow through the gap increases with 
increasing secondary mass for ijr > 0.1. For the range of mass ra¬ 
tios that we are interested in here (^ ~ 0.1 — 0.3) the leakage from 
the outer to the inner disc might be as high as 40 percent of the 
unperturbed accretion rate, further enhan cing the fossil disc mass 
dP’Orazio, Haiman & MacFadven|[2013h . Note, however, that in 
most cases these simulations, for numerical convenience, adopt 
very large values for the disc aspect ratio (H/R ~ 0.1), much larger 
than the values obtained here (H/R ~ 0.001). 

Recent 3D magneto-hvdrodvnamic (MHD) simulations by Shi 
et al. ( l2012li sueeest another mechanism that might favor an in¬ 
crease in the inner disc mass through the gap. Indeed, they show 
that in a highly ionized gas disc, where the magneto-rotational in¬ 
stability (MRI) can effectively induce angular momentum transport 
through MHD turbulence, the gas flow across the gap results en¬ 
hanced by more than one order of magnitude with respect to calcu¬ 
lations that neglect the coupling between the gas and the magnetic 
fields. 

Finally, another limitation of the model comes from the fact 
that the formulation of the problem as adopted here is strictly valid 
only if the mass ratio g <g 1. For larger mass ratios (possibly al¬ 
ready for the q = Q.\ case discussed here) both the primary and 
the secondary will rotate around the common center of mass with 
an angular velocity that takes into account both the primary and 
secondary mass. We have kept this si mple formali s m for an easy 
comparison with the previous work bv ichang et al.l ( l2010li . but we 
note here that the detailed results might be somewhat modified in 
2D or 3D, especially for the largest values of q. 


In this work we have revisited the problem of estimating the fossil 
disc mass just prior to a SMBH merger in order to explain the dis¬ 
crepancy a mong the contradictory results obta ined by different au- 
thors. While lAmitage & Nataraiaij ( l2002h and lLodato et alj ( l2009l) 
predict large values of the inner disc mass at decoupling (that thus 
result in powerful super-Eddington flaresl. lchang et al.l ( l201(]|) find 
values for the inner disc mass smaller by several orders of mag¬ 
nitude (that therefore produce significantly less powerful flares). 
These previous works consider very different scenarios for the disc¬ 
binary system and a direct comparison of the estimates of inner disc 
mass is not trivial. 

In order to find the origin of the discrepancy, we have imple¬ 
mented a ID hydrodynamic code that solves the coupled evolution 
of the disc and the binary in a self-consistent way. With our code 
we have been able to reproduce the results of Armitage & Natara- 
jan (I 2 OO 2 I) and Lodato et alj l l2009li with an extremely good agree¬ 
ment. Then, we executed the simulations performed bv IChang et alj 
(l2010l) adopting their same initial configurations and parameters: 
we found inner disc masses larger by up t o thre e orders of mag¬ 
nitude than those obtained bv IChang et al.l ( |201(]|) . As discussed in 
Section|4] we found that the origin of the discrepancy lies in the fact 
that for mass ratios larger than 0.1 the usual analytic formula for the 
tidal torque provides an unphysical contribution at large distances 
from the satellite, resulting in an artificially wide gap, and there¬ 
fore needs to be smoothed out. In our simulations, we smooth the 
torque in such a w ay to reproduce the correct gap sizes as found in 
3D simulations b y Artvmowicz & Lubovti ( Il994li . Differently from 
us, IChang et al.l 1 120101) do not perform the torque smoothing (pri¬ 
vate communication with R Chang) and thus obtain systematically 
reduced fossil disc masses and peak luminosity. 

It is worth remembering that these ID computations neglect 
gas streams that depart off the plane of the disc and that would re¬ 
quire a more realistic 3D treatment. These effects might enhance 
or decrease the actual luminosity during the merger. On the one 
hand, for a given mass qfthefossil_disCjJhe_anal^six_o£_Baruteau, 
Ramirez-Ruiz & Masset ^012l) shows that up to 80 percent of the 
inner disc mass might leak out of the secondary orbit and escape ac¬ 
cretion during the gravitational wave dominated inspiral. However, 
the limitations of such 2D computations and the unclear depen¬ 
dence of this kind of leakage on the disc-binary parameters require 
further investigations. On the other hand, leakage from the outer to 
the inner disc during the earlier disc-dominated phase would have 
the opposite effect of increasing the fossil disc mas s above what es- 
timated here. Indeed, it has been shown that MRI jShi et al 
and binary eccentricity j Artvmowicz & Lubowll 19961) are effective 
in refilling the inner disc through the gap, possibly ending with a 
further enhancement of the accretion rate on the primary BH. 

As a conclusion, taking our estimates as a lower limit to the 
fossil disc mass, even if. as obtained by Baruteau. Ramirez-Ruiz & 
Masset‘ l2012[) the~inner disc is depleted by around 80 percent, we 
would still have a super-Eddington flare in most cases, confirming 
that SMBH mergers are particularly promising for observational 
purposes. 
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